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SHOCK-TUBE FLOW ANALYSIS WITH A DIMENSIONLESS 


VELOCITY NUMBER 
By N. Muniswamy Reddy* 
Ames Research Center 


SUMMARY 


The classical shock-tube equation, for a constant-area shock tube, is 
rederived in terms of a dimensionless velocity. The equation reduces to a 
universal form so that a single graphical plot gives the solution of the shock- 
tube equation for all combinations of pressures and temperatures in the driver. 
Real-gas effects behind the shock wave are included in the solution, but the 
driver gas is assumed to be perfect with a constant ratio of specific heats. 
Specific solutions for perfect and real driven gases are discussed. All the 
thermodynamic quantities behind the unsteady rarefaction wave, ratios across 
the contact surface, wave diagram parameters and testing time in a shock tube 
are expressed in terms of the dimensionless velocity. Twin tailoring constants 
are obtained from the universal shock-tube equation, one giving the tailoring 
Mach number and the other giving the loading pressure ratio and speed -of -sound 
ratio (in the driver-driven tubes) required to achieve the tailoring Mach num- 
ber. Conditions for matching the thermodynamic quantities across the contact 
surface are given, in terms of the dimensionless velocity, and their signifi- 
cance in shock -tube performance is discussed. 


T = 


Q4* 

L-4 



Figure 1.- Wave system in an ideal shock tube. 


IMRODUCTION 


The problem of predicting the 
strength of the shock wave produced by 
instantaneously opening the diaphragm 
that separates a high and low pressure 
region has been considered extensively 
in reference 1. By using the condi- 
tions of equal pressure and equal 
velocity across the contact surface 
(C in fig. 1), one can derive an 
equation for the pressure ratio (P 21 ) 
across the moving shock in terms of 
known quantities 1 P41, y x , 


*This report was prepared while the author was a National Academy of 
Sciences' Resident Research Associate at the Ames Research Center. Dr. Reddy 
is now at Aerospace Research Laboratories, Ohio State University Research 
Foundation, Columbus, Ohio. 

1 The notation used here is the same as in reference 1; definitions of 
all parameters are given in the Notation in the appendix. Regions (l), ( 2 ), 
(3) , and (k) are identified in figure 1. 



and y 4 . For example, the equation derived in reference 1 , for a perfect 
driver and driven gas combination, in a constant-area shock tube, is 

^4 7 4 1P4&14 (fisi ~ l) 

+ [p^a/21 + 1)] 1/2 

To evaluate P 21 from equation (l) a graphical plot is necessary , since P 2 i 
is an implicit function of known quantities. Therefore, many graphical plots 
are given in reference 1 for various combinations of speed of sound and 
specific heat ratio. However, if a solution for a specific combination of 
a 4 1, P 4 1 , and 741 is desired, a new graphical plot must he generated. This 
problem becomes aggravated in combustion driven and electric arc heated shock 
tubes because different effective values of a 4 , P 4 , and 74 are obtained for 
each different run condition. Furthermore, consideration of real-gas effects 
(ref. 2) on the performance of the shock tube requires a separate graphical 
interpolation for each case. In addition, all other thermodynamic quantities 
in region (3) depend not only on the shock Mach number but also on a 4 and 
7 4 . Hence, a large series of curves for thermodynamic quantities in this 
region is needed (ref. l) . 

Hall and Russo (refs. 3 and k) recognized that equation (l) could be 
expressed in terms of two normalized parameters such that the functional rela- 
tion between them is weakly dependent on just the quantity Their work 

was limited to an ideal gas in both the driver and driven tubes. It is the 
purpose of this report to show that the shock-tube equation can be expressed 
in terms of a slightly different set of dimensionless numbers that permits one 
to include the effects of real driven gas, and moreover that permits all the 
flow quantities in different regions in the shock tube to be reduced to rela- 
tively simple functions of one dimensionless velocity number. In a sense, 
then, the work of Hall and Russo is generalized and extended here, although 
the present report pertains only to an area ratio of unity between driver and 
driven chambers, while theirs includes the effect of different area ratios in 
an ideal-gas analysis. 

In the present analysis, we rederive the shock-tube equation In such a 
form that a single graphical plot gives the solution for all combinations of 
P41 and a 4 i- Furthermore, no assumption as to the thermodynamic state of the 
gas behind the shock wave (region (2)) is made in deriving this universal 
shock -tube equation, so that it also gives the solution for a real gas behind 
the shock wave. Specific ways of obtaining the shock Mach numbers for a per- 
fect and real gas are discussed. In each case, the driver gas is treated as 
a perfect gas. However, an effective y 4 can be chosen to account for some 
of the real-gas effects in the driver. 

The following aspects of shock-tube flow are also discussed in terms of 
the new dimensionless velocity number: 

(l) Functional dependence of all thermodynamic quantities in regions (2) 
and (3), and of the wave diagram parameters. 


(P14P21) 


1 = 0 


(1) 


2 



( 2 ) Testing time. 

(3) Tailored operation. 2 The universal shock-tube equation gives twin 
universal tailoring constants. 

(k) Significance of matching the contact surface temperature. 


THE SHOCK -TUBE EQUATION 


In a shock tube, the unsteady rarefaction wave that travels into the 
driver converts the pressure and thermal energy into kinetic energy, and 
large particle velocities result behind the rarefaction wave (region (3)). 

The normal shock wave traveling into the driven gas also creates large parti- 
cle velocities in region (2) and these must match the velocity at the contact 
surface with region (3) • A general expression for this velocity, for perfect 
driver gas, is the well-known relation (ref. 1) 


u 2 = u 3 = 


2a4 


74 - 1 


[1 - (P 2 4) p4 ] 


(2) 


where the pressure, as well as the velocity, has been matched at the contact 
surface (p 2 = P3) • The maximum velocity (usually called the escape velocity) 
that can be achieved in region (3) is simply the limit of equation (2) when 


P 24 -* 0 


2a 4 

We shall define a dimensionless velocity 


( 3 ) 


u 2 u 3 

= X— = 7T- 

n U3 U3 


( 4 ) 


and in terms of this number, equation (2) becomes 


P4 B 

P 24 s (P14P21) 4 = 1 - K n 


( 5 ) 


Now, according to conservation of mass and momentum across the shock wave, 
the pressure and velocity can be expressed in terms of the shock Mach num- 
ber M s and the density ratio e = p i/ p 2 


tailoring concept is defined elsewhere in this report. 
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P 2 = Pi[l + 7^1 - e)M s 2 ] 


(6) 


= (1 - e)U s = (1 - e) ai M s 


Equation (6) may be put in the form 

*21 * ^“s 2 (l 


e 



and from equations ( 7 ) ^ (3),? and (4) 

2a 4 iR n 

M s = (1 - e)(r 4 - 1 ) 


(7) 

( 8 ) 


(9) 


or 

U s /a 4 = 2R n /[(l - e)(r 4 - l)] 
Then the shock -tube equation can be written 



which when expanded in powers of quantities less than unity becomes 



If quantities the order of |3 4 e and P 4 /y ! M s 2 are neglected compared with 

unity, which is a good approximation for Mach numbers greater than about k, 
the shock -tube equation becomes 


s n = 


(1 -E„) l/2 ^ 


v n 


(12) 


4 



where S n is the dimensionless quantity (7 X /7 4 P 4 ) (a 41 /j 7 1 P 41 ) , which we 

shall call the shock-tube number. When R n -+ 0, S n -» 00 , or ; in other words , 
a 4l -> co. When R n -* 1, S n - 0, or P 41 -* 00 . Thus, all possible shock Mach 
numbers are represented by R n between 0 and 1. The approximation neglecting 
terms of order p 4 e, that was used in deriving equation (l2) becomes invalid 
at very low Mach numbers; but normally these are not the Mach numbers of 



Figure 2.- Variation of S n vs. R n . 


interest. A plot of S n vs. R n is 
necessary to obtain R n for any given 
conditions, but a single universal 
plot (fig. 2) is sufficient to solve 
the shock-tube equation for a given 

V 

In the derivation of equa- 
tion (12) no assumption has been made 
about the thermodynamic state of the 
driven gas behind the shock wave 
(region (2)), but the driver gas has 
been assumed to be perfect. Hence, 
the value of Rn obtained from equa- 
tion (12) is unique in the sense that 
it is the same irrespective of the 
thermodynamic state of the gas behind 
the shock wave. 


SOLUTIONS OF R n 


Limiting Cases 

Perfect driven gas . - Invoking 
a 2 2 a 71 P 2 /P 2 using conservation 

of epergy across the shock wave, one 
obtains 


1 - e 


2 

7 X + 1 




Then, the relation between R n and M s (eq. (9)) becomes 


or 



Sa 4 iR n 


M s 2 » 1 


(13) 


( 14 ) 


M s ~ 5a 4 iR n 


for 




Hence, the problem of finding M s or T s from the shock -tube equation is 
solved. 

Driven gas with infinite heat capacity, e -» 0(j £ ■* l) . - In this limit 


M s = 


Real Driven Gas 

The solution for M s in the case of a real driven gas lies between the 
preceding two limits. The change in M s between the two limits is 

AM S = — (M s *) (16) 

OG 

1 

where (M g *) corresponds to the perfect gas case. Thus, the change in M s 
due to real-gas effects is typically the -order of 15 percent or less. The 
relatively small size of the effect is due to the fact that the pressure and 
particle velocity behind a shock wave are not significantly affected by real- 
gas effects. The relatively strong effects are limited to temperature and 
density. 

Exact solutions for M s can be obtained by making use of existing 
charts of real-gas normal shock properties (such as given in ref. 5 ) , along 
with the solution for R n obtained from our universal curve (fig. 2 ) . The 
value of R n determines U 2 from equation (4) and the value of M s can be 
determined from a plot of U 2 vs. M s for a given Pi. 


THE SHOCK -TUBE NUMBER S n 


2a 


’41 


- 1 


(15) 


The parameter S n in equation (12) is a function of the two important 
parameters a 4i and P 41 in a shock tube. Note that the diaphragm pressure 
ratio P 4i has far less effect on the shock-tube number than the speed-of- 
sound ratio a 41 , since S n varies as the square root of the pressure ratio 
but as the first power of the sound-speed ratio. This point has been dis- 
cussed in reference 1, but it becomes immediately evident with the present 
universal form of the shock -tube equation (eq. (12)). It can be shown that 
S n is a ratio of two limiting shock Mach numbers. 


S n = 


(M s ) 

TmJ 


P4 x~*oo 


&4 x-*x> 


(17) 
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where 


( M s)p _ 

J* 4 1 H 


2a 4l 

l7 4 " 1)11 “ e) 


(M s ) 


a 4 


N/ ^1*41 

7 1 (1 " e) 


In addition, so long as the gases in both regions (l) and (4) obey the per- 
fect gas law, it follows from the definition of S n that 


S 


n ~ 




( 16 ) 


Thus, the fundamental parameters in the analysis of shock -tube flow are the 
density ratio across the diaphragm p (eq_. (l8)) and the ratio of shock 

velocity to speed of sound in the driver gas U s /a 4 (eq. (9))- These are the 
"basic quantities represented "by our dimensionless numbers S n and Rn, and 
these basically determine the complete shock-tube performance. 


It is interesting to compare these results with those of reference 3* 

The number S n is essentially the same as Hall and Russo’s r 14 , the relation 
between them being 


S 


n _ 


2 

7 4 “ 


1 



( 18 a) 


In the ideal gas this is simply proportional to the square root of the initial 
density ratio across the diaphragm, as Hall and Russo have pointed out The 
dimension less velocity R n is somewhat different from Hall and Russo 1 s param- 
eter M s /iJP 4 i, however, even for an ideal gas: 


R 


n 


M s ideal 2 fz^±L ( \ 

a 4l gas 7 x + lJ 7 4 


( 18 b) 


It is this difference that permits us to extend the analysis to nonideal 
driven gases and to derive simple expressions for all the shock-tube proper- 
ties in terms of R n . Thus, it appears that the parameter R n is more useful 
in the shock-tube flow analysis than other dimensionless combinations. 


ACCURACY OF SHOCK -TUBE EQUATION (EQ. (12)) 


In the derivation of equation (l2) , terms of the order P^e and smaller 
were neglected. The error between the exact analysis (eq. (l)) for a perfect 
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Figure 3 . - Variation of percentage error with 
shock Mach number. 



gas and the universal solution 
(eq.. (12)) was computed at different 
shock Mach numbers and is shown in' 
figure 3. The error is about 4-. 5 per- 
cent at M s = 3 and reduces to 0.5 
and 0.25 percent at M s = 8 and 15, 
respectively, regardless of the value 
of p 4 . The magnitude of error will 
be still less for a real gas because 
e is smaller in this case than for 
the ideal gas. 


SHOCK-TUBE FLOW QUANTITIES WITH R n 


Behind the Unsteady Rarefaction Wave 

For a perfect driver gas 
(driven gas state arbitrary) the fol- 
lowing thermodynamic quantities in 
region (3) can be expressed in terms 
of R n . 3 


a 34 - (l “ Rn) 

(19) 

T 34 = (1 - Rn) 2 

(20) 

1_2 I 3 4 


= (1 - R n ) P4 

_ (21) 

1 


P 2 4 = (1 - R n ) P4 

(22) 


T 34 and P34 are plotted in figure 4. 
It can be seen that the effect of y 4 
is not very important as far as these 
quantities are concerned. 


Across the Contact Surface 

The ratios of thermodynamic quan- 
tities across the contact surface for 


3 The basic expressions for these quantities were taken from reference 1. 
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perfect driver and driven gases with M s 2 » 1 are given below with only R r 
as the variable: 


E32 - P, 


1 - Rr 


4\ R r 


T 32 = (i. 


1 - M / P 4 a i 


4l\ R,. 


32 W - Rn / \ P 4 a i 


1 ■ H n / V4 a l 


a 3 2 - 


.4 .6 

Rn 



1 ” Rn\ /^4iP4 a 


Figure 5.- Variation of B 32 with R n . 


The internal energy ratio E 32 , which 
is important in tailoring is plotted 
in figure 5- In this case, the value 
of 7 4 has a small hut noticeable 
effect . 


Shock-Tube Wave Diagram Parameters 

Speed of the tail of rarefaction 
wave . - For a perfect driver gas and 
with arbitrary driven gas thermodynam- 
ics, the velocity of the tail of the 
rarefaction wave is given as 


(X4Rv-» — 1 


If R n = l/a ,4 the velocity of the tail 
of the unsteady rarefaction wave 
becomes zero and hence the rarefaction 
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wave remains inside the driver for all times. The corresponding shock -tube 
number is 


27 ! a 4 1 



1 




(29) 


Equation (29) gives the relation between a 4 and P4 to keep the rarefaction 
wave always inside the driver. 

Tq, , X n , t s , and X s . - Using the basic equations in reference 1, we 
express the quantities r c , X c , r s , and X s (which are defined in fig. 1) in 
terms of R n : 

Perfect driver gas assumption 


. ,-0,4/2 

t c = 2(l - R n ) (30) 


X c = (a 4 - l)r c R n (31) 

Perfect driver and driven gases 
(M s 2 » 1) 


P , -a- 4 / 2 

T s = 1 _ pi/2 ^ " R n) (32) 


X s = 5R n r s (33) 


These quantities are plotted in fig- 
ures 6, 7, and 8, for various combina- 
tions of 7 4 and 7^. 
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Testing Time in a Shock Tube 

Testing time At in a shock tube is defined as the time elapsed between 
the arrival of the shock wave and contact surface at a given station along 
the driven tube. Using the relation between U2 and U s (eq. (7)) * one can 
express At in terms of R n as 

( , 4) 

x t SR n 

where x^ is the testing station along the shock tube. Once R n is obtained 
from the universal plot (fig. 2 ), At is known. It should be noted that At 
becomes very small for a real gas since e -* 0 , especially at high shock Mach 
numbers. The location for maximum testing time from figure 1 should be equal 
to X c . Then, using equations (9) and (3!)* we have 




TAILORED OPERATION IN A SHOCK TUBE 


If the reflected wave after the Interaction of a shock wave with a con- 
tact surface is a Mach wave, then the contact surface is tailored (refs. 1 
and 6) . For tailored operation, a condition of the type given "below can 
easily be derived for the perfect driver and driven gas case (refs. 1 and 6). 


and 


2 + ( 7l - l)M s g 
(37 x - DM S 2 - 2(7 x - 1) 


(36) 

(37) 


for 


Ms 2 » 1 > (^32)^ ~ 0041 

According to equation (23) for E32, the tailoring constant is given by 


(RnL = 


(aiP .) 


1/ 2 


T ~.i/2 


04 1/2 + (oiPj 1 ^ 


In terms of M s , 


(M s ) t - Sa 4 i 


(oi3.) 1/2 


a y 2 + 


In addition, the shock -tube constant for tailoring is (from eq. (l2)) 


(^n)qi - 


27^41 


(a<4i /P 4 ) 


1 

43 4 


( 7 4 - 1 ) JTjTi 


1-23 ^ 

[1 + (o4i/3 4 ) l/s ] 


C38) 


(39) 


(40) 


The twin tailoring conditions (eqs. (38) and (40)) can be represented by a 
unique point on the universal plot of the shock -tube equation as shown in fig- 
ure 2. The basic condition for tailoring is given by S n . This gives an 
explicit relation between a 41 and P 41 for tailoring which has been made pos- 
sible by the universal form of the shock -tube equation. 
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CONTACT SURFACE MATCHING 


In the tailored operation the contact surface is not totally matched in 
the sense that neither internal energy nor any of the other thermodynamic 
quantities (except P and u) are equated across the contact surface. 

(l) For a tailored operation (E 32 004 1 ) 

^'^(^TTT) Ol) 

If acoustic impedance (pa) , density, and temperature have to he matched across 
the contact surface, the following conditions in terms of temperature ratio 
must he satisfied. 

(2) For acoustic impedance matching (pa 32 = l) 

T 32 = 7 41 H 41 W 

(3) For density matching (P 32 = l) 

T32 = P 41 (^ 3 ) 

(4) For temperature matching (T 32 = l) 

T 32 =1 (44) 


It is interesting that 


(Ts 2) T ~ ( T 32) pa32=1 ~ ^32) Ps£=1 


Thus, in a tailored operation, density and acoustic impedance across the con- 
tact surface are nearly matched hut not the temperature ratio . 


For contact surface temperature matching (T 3 = T 2 ) the dimensionless 
number R n should he, according to equation (24), 


( R ri) 


TM 


(0,134) 1/2 + (8(114' 


(45) 
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The corresponding shock -tube number is 




1 


(S n ) 


n 'TM 


[1 + (pi 4 8/ai3 4 ) 1/2 3 


1 - 2(3 

4 

2B 

l/2l 4 


(46) 


It should be noted that R n and S n for tailoring (eqs. ( 38 ) and (40)) do not 
depend on the molecular weights , but R n and S n for contact surface 
temperature matching do. 

For simultaneous tailoring and temperature matching 


(%l)^ 


( R n)iM 


hence , 


P14 


^4 + ^ 

7 + 1 ~ 6 


(47) 


Then (S n )^, becomes equal to (S n )rp]y[. 

Even for different 7 combinations 0,4/8 is not very different from 
unity; hence (j, 14 (molecular weight ratio) should be nearly equal to unity 
for simultaneous tailoring and temperature matching. 

It is well recognized that in a shock-tube operation, contact surface 
mixing significantly reduces the testing time, particularly at high shock 
Mach numbers. This mixing is probably due to unequal temperature and density 
across the contact surface. Thus, it is worthwhile to match the thermody- 
namic quantities across the contact surface. This requires, according to the 
preceding discussion, that driver and driven gas molecular weights be nearly 
equal. 


CONCLUDING REMARKS 


The problem considered is that of predicting the strength of a shock 
wave produced when the diaphragm in a shock tube is opened. Without making 
any assumptions as to the thermodynamic state of gas behind the shock wave, 
but assuming that the driver gas behaves as perfect gas, we rederive the 
shock -tube equation in terms of two new dimensionless numbers. The equation 
reduces to a universal form so that a single graphical plot suffices for all 
combinations of initial conditions. The fundamental parameters that influ- 
ence the performance of a shock tube are the density ratio across the 



diaphragm, the shock velocity , and the speed of sound in the driver gas. 

All the flow parameters in a shock tube can he expressed in terms of these 
dimensionless numbers for all possible shock Mach numbers. The present 
universal shock -tube equation provides twin universal tailoring constants 
and puts the tailoring operation of a shock tube in clear perspective. It is 
hoped that the present analysis will enhance the understanding of the role of 
many variables in a shock -tube flow. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 9^035, Feb. 3, 1969 
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APPENDIX 


NOTATION 


a 

E 

M s 

P 

R n 

S n 

U s 

u 

u r 

P 

P 

At 

7 

e 

a 

0 

6 


♦ 


Jk 


speed of sound 
internal energy 

U s 

shock Mach number, — 

a x 


pressure 

dimensionless velocity (defined in eq_. (12)) 
dimensionless shock-tube number (eq.. (12)) 
shock velocity 
particle velocity 

speed of the tail of rarefaction wave 
density 

molecular weight 
testing time in shock tube 
ratio of specific heats 


P2 


(density ratio across the shock wave) 


y + 1 
7 - 1 

7 - 1 
2 / 

7i + 1 
7 A ' 1 

A (unle ss otherwise defined) 


Regions (l) , (2), (3 ), and (4) are identified in figure 1. 

Subscripts 1, 2, 3> and 4 refer to regions (l), (2 ), (3 ), and (4) as 
identified in figure 1. 
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